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Abstract 

We analysis the radial equations for massive Dirac fields in Schwarzschild black hole spacetimes. 
Different approximation formulae under the WKB scheme are developed for the transmission prob- 
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I. INTRODUCTION 



One of the most useful and efficient ways to study the properties of black holes is by 
scattering matter waves off them From the more practical point of view, the study 
of wave scattering in black hole spacetimes is crucial to the understanding of the signals 
expected to be received by the new generation of gravitational-wave detectors in the near 
future 3]. Since the linear perturbations of black holes are represented by fields of integral 
spins, the study of the scattering of wave fields are concentrated on these cases while that 
of the Dirac fields are thus less common, especially for the massive ones Recently, 
Finster and collaborators ^| have renewed the interests in that of the massive Dirac fields 
by investigating in details their evolutions in various black hole spacetimes. Here we would 
like to work in this direction by considering the scattering solutions of the radial equations 
of massive Dirac fields in spherically symmetric Schwarzschild black hole spacetimes. 

There are both numerical and analytical methods in solving the various wave equations 
in black hole scattering [lj. In this work we shall use the semi-analytic WKB approximation 
HQ, which has been proven to he very osefu, and accurate in many cases ,ike, for e X a„rple, 
the evaluation of the quasinormal mode frequencies p, |8(. For the radial Dirac equations 
we consider here, the effective potential can change from a barrier to a step and vice versa 
when the mass m or the angular momentum parameter k are varied. The WKB scheme can 
be accommodated in various ways to consider all these situations with different values of the 
energy E of the field, whether it is above or below the maximum value V m of the potential. 

In J9J], Mukhopadhyay and Chakrabarti have studied in details the radial equations of a 
massive Dirac field in the Schwarzschild black hole spacetime using a modified WKB (they 
called it the "instantaneous WKB") approximation. However, they only looked at the case 
with m = M/2, where M is the mass of the black hole, k = 1, and E m m. In this paper 
we use the standard WKB scheme instead, and we shall extend the consideration to cases 
with different values of m and k. In addition, we shall consider all possible values of E with 
E > m. In this manner we can discuss the variations of the transmission probabilities with 
respect to these various parameters. 

In the next section, we consider the reduction of the massive Dirac equation in 
Schwarzschild spacetimes into a set of Schrodinger-like equations. We also discuss briefly 
the properties of the corresponding effective potential In Section III, we look at the 
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three different WKB schemes for the cases with E 2 3> V m , E 2 m V m , and E 2 <C V m . Explicit 
approximated formulae for the transmission probability are given. In Section IV, we apply 
the formulation to calculate the transmission probabilities of the massive radial Dirac equa- 
tions for various values of E, m, and k. Conclusions and discussions are given in Section 
V. 

II. DIRAC EQUATION IN THE SCHWARZSCHILD SPACETIME 

In this section we discuss briefly the massive Dirac equation in the Schwarzschild black 
hole spacetime, including its reduction to a set of Schrodinger-like radial equations and the 
properties of the corresponding effective potentials. 

A. Radial equations 

In the Schwarzschild spacetime, 

ds 2 = -^-dt 2 + T ^dr 2 + r 2 d9 2 + r 2 sm 2 6 d<\> 2 , (1) 
r A 

where A = r(r — 2M) and M is the mass of the black hole. Consider the Dirac equation in 



this background spacetime 



Hi, 



[ 7 a e/(^ + r M )+m]^ = 0, (2) 



where m is the mass of the Dirac field, and e a M is the inverse of the vierbein e^ a defined by 
the metric g^ u , 

(3) 

with i] a b = diag(— 1, 1, 1, 1) being the Minkowski metric. 7° are the Dirac matrices 
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7°= , 7*= . ,< = 1,2,3, (4) 

I i I y id 1 J 

where o l are the Pauli matrices. Y ^ is the spin connection given by 

r M = ^[7 a ,7le/e fe ^, (5) 

where e&„ ;/Lt = d^thv — ^ u e ba is the covariant derivative of with r" being the Christoffel 
symbols. 



Here it is convenient to choose the vierbein 

/ Al/ 



(6) 



A x / 2 /r 

rsin^cos^/A 1 / 2 rsin^ shi(/>/A 1//2 rcos#/A 1//2 
rcosOcoscj) rcos^sin^ — rsin# 

— rsin#sin0 rsin#cos(/> 
Since the spacetime is spherically symmetric, one can, after some algebra, simplify the Dirac 
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equation to 



11 



12 



7°r 7A 1 / 4 d wo ma 7 - - 

7 1 7 / r 1/2 A 1/4 ^)-^-(S-L + l)^ + m^ = 0, 



+ 



A 1 / 2 <9t r 3 / 2 <9r 
where 7 is defined as 

7 = 7 1 sin# cos0 + 7 2 sin# sin0 + 7 3 cos6'. 

Also 

' B 
a 



(7) 



(8) 



(9) 



and L = f x p are the ordinary angular momentum operators. The wavefunction \1/ can be 
separated into its radial and angular parts by writing 

1 



r l/2 A l/4 



where jl^ 

with the angular parts of the wavefunction 

/ l+l/2+m" v m-l/2 
,(+) _ I V 2Z+1 r ' 



V 



(10) 



(11) 



for j = I + 1/2, and 



i+l/2-m ym+l/2 
2Z+1 ' 



l+l/2-mym-l/2 
2/+1 



a+l/2+m ym+l/2 
2Z+1 ' 

for j — I — 1/2. 

Then the radial equations 

HQ 

can be written as 
/ d 2 x 



(12) 



(13) 



dx 2 
dx 2 



+ V 1 JF = E 2 F, 
+ V 2 ) G = E 2 G, 



(14) 
(15) 



where 



with 



v lfl = ±^r + w\ (16) 

w AVy + rgg 

r 2 (K 2 + mV) + mKA/2E' 1 ' 

where 

x = r + 2 ,nn(— - l) + _ ta „- (-) . (18) 

Here we have combined the (+) and the (— ) cases, with k going over all positive and negative 
integers. Positive integers represent the (+) cases with 

k = j H — and j = l-\ — , (19) 



and 



where 




^sin(0/2) cos(#/2) 
v cos(#/2) -sin(0/2) 





(20) 



(21) 



and 3 = 1-^ (22) 



(23) 




cos(fl/2) -sin(0/2) 
sin(0/2) cos(#/2) 

From the Schrodinger-like equations in Eqs. (|14j) and (|15|) . we shall consider the scattering 
of the massive Dirac fields. Note that V\ and V2, which are related as shown in Eq. ()16|) . are 
super symmetric partners derived from the same superpotential W j^|. It has been shown 
that potentials related in this way possess the same spectra, discrete as well as continuous 
jlfil ]. Physically this just indicates that Dirac particles and antiparticles will scatter in the 
same manner around the Schwarzshild black hole. We shall therefore concentrate just on 
Eq. (|14|) with potential V\ in the following sections. 
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FIG. 1: Variation of the effective potential with the mass m (in units of M) of the Dirac field for 
k=1 and E = M. 

B. Properties of the effective potential 

Here we discuss briefly the dependence of the effective potential 



V(r, k, m, E) 

(r 2 (ft 2 + m 2 r 2 ) + mK>A/2E) 2 
A 3 / 2 ( K 2 + m 2 r 2 ) 5 / 2 



A 1/2 (k 2 + m 2 r 2 ) 3/2 + ((r - 1)(k 2 + mV) + 3m 2 rA) 
2t(k 2 + m 2 r 2 ) + 2mV + m/c(r - 1) / e] . (24) 



(r 2 (K 2 + m 2 r 2 ) + mKA/2E) 3 
on the parameters m, k, and £/. Since we shall only work with V\, but not V2, we have 
dropped the subscript of V. 

First, its dependence on m is showed in Fig. ^ with energy E = M and with k — 1. Note 
that we have shown V as a function of x in the figure. For small values of m, the potential 
is in the form of a barrier, with the asymptotic value 

V(x -> 00) = m 2 . (25) 

As m is increased, the peak of the potential also increases but does so very slowly. Eventually, 
the height of the peak is lower than the asymptotic value m 2 . When m is increased further, 
the peak disappears altogether, and the potential barrier turns effectively into a potential 
step. 

The effective potential V also depends on the energy E. However, as shown in |10], the 
general behaviors of the potential remain the same as E is increased from its minimum value 
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FIG. 2: Variation of the effective potential with k for the Dirac fields of m = 0.2M and 0AM. 

m. Indeed, from the form of the potential in Eq. (|24jl . E appears all in the denominators. 
The terms involving E can never get large enough to change the general behaviors of the 
potential since E cannot be smaller than m. 

Finally, in Fig. |21 we show the dependence of the effective potential with the angular 
momentum quantum number k. We see that as k increases, the behaviors of the potential 
approach to that of the massless one, regardless of the mass of the field. This can be seen 
by taking the limit \k\ —>■ oo in Eq. (|24|) . 

V(\k\ - oo) « (26) 

which is independent of m. 
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III. WKB APPROXIMATIONS 



In this section we consider the radial equations obtained in the previous section as a 
general one-dimensional quantum mechanical problem, 

{-1& +V )* = **< < 27 > 

with the asymptotic values of the potential 

V(x -> oo) = m 2 and V(x -> -oo) = 0. (28) 

Suppose an incident wave comes from the right (x = oo). Then the boundary conditions for 

ip are 



^(x^oo) = e - WE2 - m " x + Re iVE ' 2 - m ' x (29) 
i/i(x -> -oo) = Te~ iEx , (30) 

where R and T represent the reflection and transmission coefficients, respectively. 

Since exact solutions for Eq. (|27j) are usually hard to find, we have to resort to approxi- 
mations. Here we adopt the WKB approximation, which has been proved to be extremely 
useful and sometimes to be more accurate than expected. We try to develop approximate 
expressions for the transmission probability 



for the whole range of the energy E, including the cases with E 2 ^> V m , E 2 w V m , and 
E 2 V m , where V m is the maximum value of the potential. In the case of a barrier, V m will 
be the peak value of V, while in the case of a step, V m = m 2 . 

A. E 2 > V m 

For E 2 ^> V m , the standard WKB form of the wavefunction can be given by 

tP( x ) = C+W + (x) + C7_W_(x), (32) 

for — oo < x < oo, with 

ttt I \ 1 ±i dx'p(x') /oo\ 

W±{x) = e Jx o (33) 

/p(x) 
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being the WKB wavefunctions. Here p(x) = \JE 2 — V(x), xq is some fixed reference point, 
and C + and C_ are constants to be determined from the boundary conditions. Eq. (J3*2j) is 
a useful approximation for the wavefunction as long as the validity condition 

V 



2p 3 



<1, 



(34) 



where V = dV/dx, is fulfilled. 

From the boundary condition at a; — > — oo (Eq. (|3Uj) ). we have C + = 0, while from the 
condition at x — > oo (Eq. we have 



C_ = (E 2 - m 2 ) l/i e~ WE2 ~ m2 ^og*/^^'^')-^-™ 2 ) 



(35) 



Therefore, the transmission coefficient 



E 

' E 2 — m 
E 2 



2\ 1/4 



,i(E-y/&-m*)x<i jf™ dx'(p(x')-VE-'-m^ i dx'(p(x')-E) 



(36) 



and the transmission probability in Eq. (JBTj) becomes 



T = 1. 



(37) 



This is of course consistent with the classical result in which the particle moves to the left 
without rebounding. However, quantum mechanically T is not exactly equal to 1, there is 
always a small probability for reflection. Thus to get a better approximation on T in this 
quantum mechanical situation, we need to extend the WKB scheme used here. 
Suppose we write the general solution of Eq. (j77|) as 



ip = C+{x)W+(x) + C-{x)W-(x), 



(38) 



where C + (x) and C-{x) are now functions of x. If we also take 

C' + (x)W + (x) + C'_{x)W-{x) = 0, 
and together with Eq. (|38|). we can obtain the differential equations for C±, 



C' ± {x) = T i 



' p"(x) _ 3(p'(x)y 
Ap(x) 2 8p(x) 3 



c±(x) + c T ( x y •>*<> 



=F2i r dx'ptx') 

J xn 



(39) 



(40) 



9 



From the boundary conditions in Eqs. ()29|) and (|30|) . we have 
C+(-oo) = 0, 

CJoa) = (E 2 - m 2y/4 e -i^^ J™ dx 'W E2 - v ( x ')~^ E2 - m ' 2 ) ^ 



(41) 
(42) 



Finally, with these boundary values, the differential equations in Eq. (|4f)|) can be turned into 
integral equations, 

,( P "{x') 3(p\x')r 



, ( P "(x') 3( P 'lx')r 



„ i -2i r dx" P (x") 

C-(x)e Jx o 



C-(x) = C_(cx>) — % / dx 



4p(x') 2 8p(x') 3 



C(x') + C + (x')e J *o 



n 2i L d *"p(*") 



(43) 



(44) 



To the lowest order, we just substitute the boundary values of C± into the right hand side 
of the above equations, 



C+(x) ss -iC. 
CJx) « C_(oo 



V J J-oo \Ap(x') 2 8p(x') 3 J 



\4p(x') 2 8p(x') 3 



Under this approximation, the reflection coefficient becomes 
R 

with 



{E 2 -m 2 ) l l i 



-{oo) J^ax ^ Ap{x , )2 8p{x , )3 )e 



C+(oo) = -iC. 

as given in Eq. ([45)1 . The reflection probability in this approximation becomes 



(45) 
(46) 

(47) 
(48) 



K = \R\ 



50 dx' ( p " {x ' ] - e -<-w> 

\Ap(x') 2 8p(x') 3 J 



(49) 



while the transmission coefficient is given by 



T = 1-TZ 
= 1- 



, f P "(x') s(p'(x')) 2 \ » r' 



dx 

oo ^^(a;') 2 8p(a; / ) 



e x o 



(50) 



This constitutes our WKB approximation for i? 2 3> V m as long as the validity condition in 
Eq. (|34j) is satisfied. Note that we can also obtain, in this approximation, the wavefunction 
by substituting C + (x) and C_(x) in Eqs. (J43j) and (jH3), respectively, into Eq. 
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B. E 2 « V m 



In the case of a potential barrier, when the energy E of the field is decreased to such an 
extend that E 2 is close to the peak of the potential, the validity condition (Eq. (|34j) ) will 
no longer be satisfied near x m , the position of the peak. This situation can be remedied by 
representing the part of the potential parabola, while maintaining the WKB 

solutions on either side of it. Exact solutions can be found for the parabolic potential, and 
then the approximate solution for the whole range, — oo < x < oo, can be obtained by 
matching the three solutions across the intertwining regions j^, 0| . 

For the part of the potential near x m , we can write, in the parabolic approximation, 

V(x) ^V m + ^V"(x m )(x - x m )\ (51) 

and 

E 2 -V(x) « (E 2 -V m )+X(x-x m ) 2 

VX(z 2 + £ 2 ) E 2 > V m 

^ v\(z 2 - a e 2 < v m 

where 

X = -V"(x m )/2, z = X 1 /\x-x m ), and ^=\^^H. (53) 
Note that for 

\E 2 - V m \ 1/2 

Z = ±^ X h2 = X m ± . (54) 

When E 2 is smaller than V m , X\ and X2 are just the turning points with E 2 = V{xxp)- In 
any case for E 2 w V m , we can divide x into three regions: (I) x > x±, (II) x\ > x > x%, and 
(III) x < X2- For regions (I) and (III), we still use the WKB form of the wavefunction, 

tfjj(x) = A + W + (x) + A^W^(x), 
iPui{x) = BW + (x). (55) 



For region II, we have, in the parabolic approximation, the Schrodinger equation in Eq. 
can be written as 

J 2 ^ 

(z 2 ±e)^ = 0, (56) 



dz 2 

with the general solution 



ijjj(z) = aD , ^(V^/V) ^(-V2e^z), (57) 

2+ 2 2+ 2 
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where D u (t) are the parabolic cylinder functions. 

Taking into account of the boundary conditions in Eqs. (j22j) and ()30|) . and matching the 
wavefunctions in different regions across x = xi and Xo, we can solve for the constants A+, 

nn 

A_, a, /?, and B |6J,|8|]. One can then obtain 

T J i/d + ^>v„ 

(l/(l + e"« 2 ) £ 2 <K„ 
These two cases can be combined by writing 

T ~ i + eTr^-i? 2 )/^' ^ 
which is the WKB approximation we shall use for E 2 « K m . 

C. £ 2 < y m 

When the energy of the field is lowered further, the turning points x\ and Xi will move 
far apart in such a way that the parabolic approximation is no longer valid. In this case, 
one can divide the x-axis into five regions: (I) x > x\, (II) x ~ xi, (III) X\ > x > X2, (IV) 
x ~ X2, and (V) x < x%- For regions (I), (III), and (V), we still use the standard WKB 
wavefunctions, 

$j{x) = A + W+{x) + AJW-(x), 
^jjj(x) = B + W + {x) +5_iy_(x), 
if> v (x) = CW + {x). (60) 

For region (II), we use the linear approximation Q|, 

V(x)^V{x 1 ) + V'(x 1 ){x-x 1 ), (61) 

and 

E 2 - V(x) « t i 2 l /3 z 1 , (62) 

where 

jJi = —V'(xi), and zi = ^ 3 (x — Xi). (63) 
Then the Schrodinger equation becomes 

% + Z ^ = °' (64) 
12 



with the general solution 

iP n ( Zl ) = aAi(-zi) + (3Bi(~zi), (65) 

where Ai(t) and Bi(i) are Airy functions. One can consider region (IV) in a similar way to 
obtain 

ij IV (z 2 ) = 1 Ai(z 2 ) + 5Bi(z 2 ), (66) 

where 

A*2 = -V(x 2 ), and z 2 = n 2 (x - x 2 ). (67) 

Again matching the boundary conditions at x = ±oo and the wavefunctions across x = x\ 
and x 2 , one can obtain the constants A + , A_, B + , £>_, C, a, (3, 7, and 5. The transmission 
probability is then given by 

r = e - 2/ ;>Vn^ (68) 

This is our WKB approximation for the cases with E 2 V m . 

IV. TRANSMISSION PROBABILITIES FOR THE DIRAC FIELD 

Using the WKB approximations outlined above we shall calculate in this section the 
transmission probabilities for the Dirac field in the radial equation (Eq. (fH|) ) with the 
potential as given by Eq. (jUJ). 



A. Massless cases 

First, we consider the massless cases, with the potential 

V(r, K ) = ^^[| K |A 1 / 2 -(r-3)], (69) 

where we have used the mass M of the black hole as a unit of mass and length to simplify 
the notation so that A = r(r — 2). To proceed we first consider the case with k — 1. Then 
we have 

^ = 0.216. (70) 

Hence, for E ^> 0.216, we can use the formula in Eq. (|5()j) to evaluate the transmission 
probability. To check the validity condition, we plot |V/2p 3 | for various energy E in Fig. 01 
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FIG. 3: Validity condition for various values of E of the massless Dirac fields with n = 1. 



From it, we see that for E = 0.35, the maximum value of this quantity, 

V 



2p 3 



0.1. 



(71) 



The formula in Eq. (J50|) can therefore be a good approximation for energy E larger than 
around 0.35. The transmission probabilities T for 0.35 < E < 1 calculated from this formula 
is plotted in Fig. HJ 

For E smaller than 0.35, one can use the parabolic approximation for the transmission 
probability as given in Eq. ([59)1 . The transmission probabilities for < E < 0.35 calculated 
from this approximation is also plotted in Fig. 0] Near the base of the barrier, the parabolic 
approximation is no longer valid because the turning points are far apart. This is apparent 
from the fact that the curve tends to around 0.0377 instead of zero as E goes to zero. 

When the turning points are isolated, one can use the approximation given in Eq. ()68|) for 
the transmission probability T . The transmission probabilities for < E < 0.216 calculated 
from Eq. (JoTJj) are also plotted in Fig. 0] When E is close to \/V m = 0.216, the formula in 
Eq. (|59|) cannot be trusted as we can see that T — > 1 as E —>■ V m in this approximation. 

In Fig. HJ we see that the solid and dashed lines calculated from these two approximations 
overlap at around E = 0.154. Thus for E > 0.154, we should take the result with the 
parabolic approximation (solid line) and for E < 0.154, we should take that of the tunneling 
approximation (dashed line) instead. While in the region with E ~ 0.154, we can extrapolate 
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FIG. 4: Transmission probabilities of the massless Dirac field with k = 1 in the various WKB 
approximations for E 2 3> V m (dotted line), E 2 ~ V m (solid line), and E 2 <C V m (dashed line). 
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FIG. 5: Transmission probabilities T of the massless Dirac field with k = 1 to 5. 

smoothly between these two curves. Combining with the result for E > 0.35, we can obtain 
a curve for the transmission probabilities in the entire region, < E < 1. This is plotted in 
Fig. El 

In addition to the curve for k = 1, we also plot the transmission probabilities for k = 2, 
3, 4, and 5 in Fig. These curves are obtained in the same procedure as outlined above 
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for k = 1. They are similar to each other and shift to the right as k increases. The energy 
at which T = 1/2 occurs is \/V m which increases as the peak of the barrier gets higher and 
higher when k is increased. 

B. Massive cases 

The situations with nonzero m are more complicated because the potentials can change 
from a barrier to a step or vice versa as the parameters m and k are varied as shown in 
Figs. H and El With k = 1, the potentials are barriers for m = 0, 0.1, and 0.2. One can 
use the relevant approximations for different values of the energy E, along the same lines 
as in the massless cases in the last subsection. The results are plotted in Fig. |U1 However, 
for m = 0.3, 0.4, and 0.5, the potentials are steps. The only approximation one needs is for 
E m. We also check the validity condition and the result indicates that the approximation 
is useful all the way down to energy value very close to m. The results for these masses 
are also plotted in Fig. |H1 For the potential steps, the transmission probabilities are almost 
equal to 1, which are consistent with the classical results. 

We also see from the various diagrams in Fig. |U] that the variations of the transmission 
probabilities T with m are numerically very small. In order to see the changes in more 
details, we plotted T versus m in Fig. for E = 1 and k = 1. Since the transmission 
probabilities in these cases are very close to 1, we take the logarithmic of T in the plot. 
From this figure we see that T first decreases from m = 0, attends a minimum around 
m = 0.14, and then increases as m is further increased. 

The variations of T between m = 0.3 to 0.5 are quite unexpected. This is because 
with these values of m, the potentials are in the form of steps with larger and larger step 
height when m is increased. For simple potential steps, it is known that the transmission 
probabilities decrease as the steps get higher and higher for fixed energy. We thus see that 
for the black hole effective potentials, the transmission probabilities are not determined only 
by the height of the potential but also the overall shape of it. 

In Fig. El -£7 = 1 and is much larger than -y/Kn- For energy values closer to y/Vm, the 
variations can be quite different. As shown in Fig. |H1 for E = 0.2, that is, when E 2 is near 
the peak of the potential, T increases first as m is increased from 0, attends a maximum 
value around m = 0.09, and then decreases. The same trend is found for E = 0.1 where it is 
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FIG. 6: Transmission probabilities T of the Dirac field with k = 1 and m = to 0.5. 
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FIG. 7: Variation of the logarithmic of the transmission probabilities of the Dirac field (E = 1 
and k = 1) with m. 
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FIG. 8: Variation of the transmission probabilities of the Dirac field with m for E = 0.1 and 0.2. 

well below the peak of the potential. However, the maximum value in this case is at around 
m = 0.05. 

As we can see from the above discussions, the transmission probability T in general 
increases with the mass m when m is large enough that the effective potential is in the form 
of a step. However, when m is smaller and the effective potential is in the form of a barrier, 
the variations of T with m can be quite complicated. When E is large and well above the 
peak, T first decreases and then increases when m is increased from m — 0. When E is 
small with its value near or well below the peak of the potential, the variation is reversed, 
that is, T first increases and then decreases when m is increased from m — 0. 
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V. CONCLUSIONS AND DISCUSSIONS 



We study the radial equations of the massive Dirac field in the spherically symmetric 
Schwarzschild spacetime. Using the WKB approximations and the appropriate connection 
formulae, we are able to give semianalytic formulae for the transmission probability T of the 
radial wavefunction with the energy E 2 ^> V m , E 2 m V m , and E 2 <C V m , V m is the maximum 
value of the effective potential. For the massless cases, we find that, as shown in Fig. the 
variations of T with the energy E and the angular momentum number k are similar to that 
for the scalar case ^, ^|. Since the potentials are in the form of barriers and the heights 
of the peaks of these barriers increase with k, T for fixed E will thus decrease as expected. 
This means that waves with lower angular momenta but with fixed energy will be absorbed 
more easily by the black hole. 

The massive cases are more complicated. The effective potentials in these cases change 
from barriers to steps when the mass of the field is increased. When the potential is still 
in the form of a step and the energy E of the field is well above the maximum value of the 
potential, the transmission probability T decreases with m. Thus higher mass fields will get 
absorbed by the black hole more easily. However, when m is decreased to such an extend 
that the potentials become barriers, the variation trend changes as shown in Fig. [7| At some 
value of m (around 0.14 for E = 1 and k — 1), T turns around and increases when m is 
further decreased to 0. Therefore, we see that the variations of T with m is complicated 
when the potentials are in the form of barriers. This is also true when E is smaller with its 
value near or well below the peak of the potential as shown in Fig. |SJ 

After calculating the transmission probabilities for the radial wavefunctions, one can 
evaluate the corresponding phase shifts and cross sections in various scattering situations 
Q]. In the semiclassical limit one can also study the interesting phenomena such as black hole 
glories 1^1 . orbiting and spiraling scatterings lfj by deriving the semiclassical deflection 
function [2£|. We hope to further investigate these issues in our future publications. 

Quite peculiarly for black holes, one can also study the absorption cross sections for 
various wave fields. There has been quite a lot of interest in these cross sections in relation 
to the higher- dimensional black holes in string theory, especially the low-energy absorption 
cross sections Q. However, the WKB approximations that' we use in this work is not 
adequate at low energy, that is, when Emm. One can nevertheless improve the WKB 
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approximations in these threshold situations when E is near the top of the step or when it is 



near the base of the barrier 



22, 



This improved approximations may therefore provide an 



alternative to the usual method |3| in obtaining these low-energy absorption cross sections. 
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